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Abstract. A straightforward formalism to evaluate the impact of Type la 
Supernovae (SNIa) on large scale astrophysical issues is presented, together with 
analytical formulations for the SNIa rate following an instantaneous burst of Star 
Formation (SF), for the variety of SNIa progenitors. Some applications of the 
parametrized formalism are illustrated. The observations seem to favour the 
Double Degenerate (DD) systems as SNIa precursors. 



1. Introduction 

Type la supernovae are an important source of Iron and energy to the interstellar 
medium, and for this reason the evolution of their rate over cosmic times has 
profound implications on a variety of astrophysical issues, including the chemical 
evolution of galaxies, the Fe content of clusters of galaxies, the evolution of 
the gas flows in Ellipticals (Es). For example, the delayed release of Fe from 
SNIa with respect to the prompt release of a elements from Type II supernovae 
implies that the shorter the formation timescale of a stellar system, the higher 
the a/Fe abundance ratio recorded in its stars. Following this argument, short 
formation timescales have been inferred for E galaxies from the supersolar a 
to Fe ratios traced by the Mg and Fe indices in their spectra (e.g. Matteucci 
1994). As another example, the interplay between the rate of mass return from 
evolved stars in Es and the rate at which this matter is heated by SNIa explosions 
crucially determines the dynamical evolution of the gas in these galaxies (Renzini 
1996), i.e. whether from inflow to outflow or from outfow to inflow. 

Our ability to address these (and other) issues is hampered by the uncer- 
tainty which still affects the SNIa progenitors. Current models include Single 
Degenerate (SD) and Double Degenerate (DD) systems, both dealing with a 
close binary in which the primary component (mi) evolves into a CO white 
dwarf (WD). When the secondary (771-2) evolves off the MS, Roche Lobe Over- 
flow (RLO) occurs: in the SD model, the WD accretes and grows in mass; in the 
DD model, the envelope of the secondary is lost from the system, leaving a close 
double WD, which eventually merge due to the emission of gravitational radia- 
tion. In both cases, explosion may occur either when the accreting WD reaches 
the Chandrasekhar limit and C deflagrates (Ch), or when a He layer of ~ 0.15 
Mq has accumulated on top of the CO WD so that He detonates (S-Ch) (Hille- 
brandt &; Niemeyer 2000) . The various models correspond to markedly different 
secular evolution of the SNIa rate, with great variance of the consequences on 
related issues. 
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To assess the impact of the SNIa rate on any astrophysical problem it is 
important to (i) adequately characterize the progenitor model and (ii) correctly 
couple this model to the SF history of the specific system. In this paper I will go 
through these steps, and show how the scenarios for the SNIa progenitors can 
be constrained by their impact on large scales. To do that I will use analytical 
relations for the SNIa rate for both SD and DDs derived in Greggio (2004), as 
a convenient parametrization of the progenitors' models. 



2. The Model 

An effective formulation of the SNIa rate for the modelling of stellar systems 
rests upon the definition of two key quantitites: the distribution function of the 
delay times fi a (r), where r is the time between the birth of a star and its death 
as a SNIa, and the number fraction of SNIa events out of one stellar generation, 
Aj a . Denoting with t\ and r x the minimum and maximum delay time, and if A\ a 
is constant, the SNIa rate at epoch t in a system with a SF rate ip is: 

/■min(t,T x ) 

ni a (i) = k a ■ A la ■ / - t) ■ /i a (r) dr (1) 



where k a is the number of stars per unit mass in each stellar generation (e.g. 
k a = 2.83 for Salpeter IMF between 0.1 and 120 Mq), considered constant in 
time. 

Eq-lfl]) gives the dependence of the observed SNIa rate on the SF history: 
for systems with an almost constant vp, like late type galaxies, at late epochs 
(t~r x ): 

r&(t) ~ W • k a ■ A Ia (2) 

since /ia(i~) is normalized to one. Thus, the current SNIa rate in late type 
galaxies essentially gives information on the realization probability of the SNIa 
scenario, since it results from the contribution of systems with all possible delay 
times. By approximating {ip) with the ratio of the galaxy mass over its age, and 
assuming an M^/L-q ratio of the order of unity, Eq.(J2J) requires A\ & ~ 10~ 3 to 
match the observed value of 0.2 SNUs (Cappellaro, Evans, & Turatto 1999). 

In early type galaxies, where most SF has occurred in an initial burst of 
relatively short duration At, the late epoch SNIa rate is: 

hf a {t) = M* • k a ■ Ai a • (fu)t-At,t (3) 

i.e. it is proportional to the value of /i a at delay times close to the age of the 
galaxy, and to the total mass of stars formed in the burst (M*.). 

If Ai a and k a are the same in early and late type galaxies, the ratio between 
the current SNIa rates measured in SNU in the two galaxy types is: 

t-. rfa^SNU (M+/L B ) E (fla)t-At,t , ,x 

n Ia,SNU ( M */ L B) L (fl*)n,t 
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Using Maraston (1998) solar metallicity models, at an age of 12 Gyr a system 
formed at a constant rate has M+/L-Q ~ 2, while a single burst stellar population 
has M*/L B — 13. Since 7^-snu ~ 1 (Cappellaro et al. 1999), Eq. (gj) implies: 

\/la)t— At,t n -, r / c \ 

— — 0.15 (5) 

\Jla.)n,t 

showing that, on the average, the distribution of the delay times must be de- 
creasing with r, i.e. the majority of SNIa precursors are relatively short lived. 

The theoretical SNIa rate can be obtained with numerical simulations, 
which follow the evolution of a population of binary systems, under some pre- 
scriptions concerning the outcome of the mass transfer phases (e.g. Han 1998). 
Typically, the predicted A\ a ranges within 10~ 4 and 10~ 3 , while /i a is charater- 
ized by an early steep rise, followed by a decline at later epochs. These results 
depend on a variety of input parameters, whose role is difficult to gauge, and 
therefore they are hardly suitable for a thorough exploration of the parameter 
space. On the other hand, the shape of the /i a (T) function can be characterized 
on the basis of general considerations founded on stellar evolution. Along this 
line, Greggio (2004) derives analytical relations for the /i a function for both the 
SD and the DD model, briefly described below. 




Figure 1. The SD model: 71(7712) is derived from the distributions n(mi) oc 
m^ a and n{q) oc q 1 where q = 7712/7711, and summing over the relevant combi- 
nations. The right panel shows the ff a D function on arbitrary units, as from 
Eq. ©. r is the MS lifetime of the mass 7712. e is the ratio between the ac- 
creted mass and the envelope mass of the donor. Curves are shown for e = 1 
and e = 0.5 (labelled). For the Sub-Ch models, it is assumed that m\ > 3. 

In the SD model, the delay time is equal to the MS lifetime of the secondary, 
plus the duration of the accretion phase, which is negligible. Thus: 

fi D {T)^n(m 2 )-\m 2 \. (6) 

where \rri2\ is the rate of change of the mass whose MS lifetime is r, and 71(7712) is 
the distribution function of the secondary masses in those systems which produce 



4 



Greggio 




8 8.5 9 9.5 10 8 8.5 9 9.5 10 



Figure 2. The DD model: dependence on r njX (left), and on the slope of 
the distribution of the gravitational delays (right). See text for more details. 



SNIa events. Fig. ^ shows the two factors in Eq.© (left panel), and their 
product (right panel). The shape of the n{m,2) function reflects the limitations 
imposed on the components masses in order to give rise to a SNIa. It is assumed 
that mi varies between 2 and 8 Mq, since most CO WD come from this range. 
The lower boundary originates a relative lack of systems once mi drops below 2 
Mq. The last steep decrease, instead, reflects the need of building a WD up to 
the Ch limit: as mi decreases, so does its envelope, that is the donated mass. 
Then, at late epochs, only the most massive WD will end up as SNIa, which 
implies a lower limit of mi > 2M & . Below rri2 = 2M , 12(1712) is very sensitive 
to a, 7, e (defined in the caption), which determine the volume in the parameter 
space of SNIa progenitors. The resulting /i a function in shown in the right panel: 
notice the dramatic drop at late epochs which stems from the requirement of 
building up to the Ch mass. S-Ch exploders (dot dashed line) do not show this 
drop, the only requirement being that the CO WD is more massive than 0.7 Mq. 
If additional limitations are imposed on 1712, the /j a function is non zero only 
in the range of r equal to the MS lifetimes of the secondaries with appropriate 
mass. For example, if systems with 7712 < 2M Q evolve into a common envelope, 
/la = f° r T > 1 Gyr. On the other hand, in order to explain the occurrence 
of SNIa in Es with the SD model, low mass secondaries must be allowed for. 

In the DD model, the delay time is equal to the MS lifetime of the secondary 
(r n ), plus the gravitational delay (r gw ). Considering only CO DDs, 772-2 is (in 
most cases) heavier than 2 Mq, which implies r n < 1 Gyr; r gw instead can span 
a very wide range, being extremely sensitive to the separation A of the DD. For 
example, for a (0.7+0.7) Mq system, r gw goes from 14 Myr to 18 Gyr when A 
increases from 0.5 to 3 Rq. Since the separation of the progenitors of the CO 
DDs ranges from some tens to some hundreds Rq, a great degree of shrinkage 
must occur at the RLO phases in order to produce explosions within a Hubble 
time. /j^ D can be derived as a modification of ff®, with early explosions from 
systems with short r n AND r gw , and late explosions from systems with long t^. 
The distribution of r gw will be very sensitive to the distribution of A within the 
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range of a few Rq. Motivated by literature results, Greggio (2004) considers two 
possibilities: (i) the WIDE DDs, in which the first RLO corresponds to a modest 
shrink of the system (as, e.g., in Nelemans et al. 2001); (ii) the CLOSE DDs, 
similar to the classical recepy (as, e.g., in Han 1998) . The two channels imply 
different characterizations of the distribution of the gravitational delays. Fig. 
121 shows the resulting functions for some choice of the major parameters: 
T nx , the MS lifetime of the least massive secondary still progenitor of a SNIa; 
/3 a and f3 g which are the exponents of the two power laws assumed to represent 
the distribution of A for the WIDE DDs, and of r gw for the CLOSE DDs l . 

The fR function shows an early rise, a wide maximum and a late epoch 
decline. The width of the maximum is ~ r njX . Relative to the CLOSE DD, the 
WIDE DD scheme of evolution leads to flatter distributions, and the more 
systems with short (i.e. steeper /3 a or (5 g ), the steeper the early rise as well 
as the late epoch decline. 

3. Constraining the SNIa Progenitors and Implications 

In summary, stellar evolution arguments strongly support a distribution function 
of the delay times characterized by: 

• an early maximum, soon after the most massive SNIa progenitor explodes; 

• a relatively flat portion, up to a delay time equal to the MS lifetime of the 
least massive secondary in binaries which end up as SNIa; 

• a decline phase, which for the SD Ch models becomes very steep at late epochs, 
due to the requirement of building up the Chandrasekhar mass. 

In spite of the similar shape, the different models have very different as- 
trophysical implications: this can be appreciated by looking at the left panel in 
Fig. El where the various models are plotted normalized at 12 Gyr, so that they 
all reproduce the currently observed SNIa rate in Ellipticals. Notice that f\ & 
is proportional to the SNIa rate following an istantaneous burst of SF, not too 
far from what is expected in Es, whose stars are mostly old. The SD Ch model 
predicts a huge variation of the SNIa rate, which at early epochs would have 
been ~ 2000 times greater than now. Similarly, the S-Ch model corresponds 
to a large excursion of the SNIa rate, but, after the maximum, the rate drops 
rather quickly. Finally, for the DD models the initial peak is less strong. 

The various models are characterized by different values of the ratio in 
Eq.©, and they correspond to vastly different total Fe production, as shown 
on the right panel of Fig|31 According to Renzini (2003), the Fe Mass-to- light 
ratio in Clusters of Galaxies, which are dominated by Es, is ~ 0.015: if 2/3 
of this is provided by SNIa, then M^/Lb ~ 0.01. This, admittedly rough, 
estimate is reported in Fig|21 together with the empirical value for the ratio in 
Eq.©. It appears that, once normalized to reproduce the current rate in Es, SD 
models greatly overproduce Fe in galaxy Clusters AND overpredict the current 
SNIa rate in late type galaxies, where all the delay times are sampled. The 
data require a less dramatic secular evolution of the SNIa rate, similar to the 
behaviour of the DD models. 



1 If SNIa originate from DD systems with approximately the same mass, /3 g ~ —0.75 + 0.25/3 a 
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Figure 3. Normalized /i a (r) functions (left) and their global properties 
(right) for the various progenitors models. The total Fe production per unit 
luminosity of the parent galaxy is computed by adopting a current SNIa rate 
in Es of 0.2 SNUs, and a release of 0.7 M Q of Fe per SNIa. The numbers in 
parenthesis are respectively r njX and f3 a (CLOSE, DD-C) or (3 a (WIDE, DD- 
W). All the plotted models assume a = 2.35,7 = 1 to derive the distribution 
of the secondary masses in the primordial binaries. 



The different models in Fig. |3] imply different timescales over which, fol- 
lowing a burst of SF, 1/2 of the total SNIa explosions occur. Although the DD 
model is characterized by a less dramatic evolution of the SNIa rate, since most 
of the explosions take place within the first Gyr, supersolar a/Fe ratios still 
imply SF timescales within 1-2 Gyr. Regarding the dynamical evolution of the 
gas in Es, detailed modelling is needed, since at late epochs the SNIa rate from 
DDs and the rate of mass return decline at a similar pace (oc r~ 13 ). 

In summary, taking advantage of an analytical fomulation of the distribution 
of the delay times of SNIa progenitors, I have outlined a straighforward way to 
estimate their impact on the large scales. A more detailed elaboration of the 
issues considered here, as well as other applications, will be presented elsewhere. 

References 

Cappellaro, E., Evans, R., & Turatto M. 1999, A&A 351, 459 
Greggio, L. 2004, A&A, submitted 
Han, Z. 1998, MNRAS 296, 1019 

Hillebrandt, W., & Niemeyer, J. C. 2000, ARA&A38, 191 
Matteucci, M. F. 1994, A&A 288, 57 
Maraston, C. 1998, MNRAS 300, 872 

Nelemans,G., YungelsonRR., Portegies Zwart,S.F., & Verbunt, F. 2001, A&A 365, 491 
Renzini, A. 1996, in IAU Symp. 171, New Light on Galaxy Evolution, ed. R. Bender 
& R.L. Davies (Dordrecht:Kluwer), 131 

Renzini, A. 2003, in Clusters of Galaxies: Probes of Cosmological Structure and Galaxy 
Evolution, cd. J.S.Mulchaey, A. Dressier, & A. Oemler (CUP), 261 



